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SUMMARY 


Thin flat slabs of scintillator are a useful means of measuring 
angular location of gamma-ray fluxes of astronomical interest. A 
statistical estimate of position error has been made of two scintillator 
systems suitable for gamma-ray burst location from a balloon or 
satellite platform. A single rotating scintillator with associated 
flux monitor is compared with a pair of stationary orthogonal 
scintillators. Position error for a strong burst is of the order of a 
few arcmin if systematic errors are ignored. 



INTRODUCTION 


Slabs of scintillator may be used as directional gamma ray detectors 
for a number of studies in astronomy including measurement of cosmo- 
logical gamma ray anisotropy, location uf discrete sources, measurement 
of atmospheric anisotropy, and location of gamma ray bursts. 

We examine here only the last of these, and attempt to estimate the 
position error of scintillator systems of reasonable dimensions used to 
locate these transient phenomena. 

Five years ago the gamma ray bursts were discovered by detectors 
on V^la spacecraft. Since then some 40 bursts have been detected by 
at least 9 spacecraft in addition to several Vela satellites. The 
bursts are a unique phenomenon with no parallels in other energy 
regions or time scales of astronomy. Compared to other astronomical 
sources of gamma rays the bursts have large flux, a hard spectrum and 
duration of less than one second to some tens of seconds. They have 
not yet been identified with any optical or x-ray object. Crude position 
measurements with uncertainties of 5°-10° are available for less than 
a dozen bursts and so far indicate an apparently random distribution. 

The next major step in the understanding of these extremely 
energetic processes requires their position location to a degree of 
accuracy which allows identification with an object previously observed 
in some other spectral region. Candidate instruments and techniques 
which have been proposed for this purpose, include long baseline timing,-** 

r\ ^ 

the x-ray shadowgraph, an array of one dimensional Dicke cameras, and 



an active anti-collimator.^ We examine here an alternative experimental 
arrangement using the directional sensitivity exhibited by a thin flat 
scintillator to a parallel beam of gamma rays. The general case of 
response of such shaped scintillators has been treated in detail by a 
number of authors, most recently Trombka et. al.~* (1975). We have 
selected here a scintillator of dimensions suited to the purpose of 
locating the direction of the gamma ray bursts from a small satellite. 

We have examined two configurations in detail. The first consists 
of a planar scintillator rotating rapidly about an axis in its own 
plane. A single rotating scintillator would locate gamma ray bursts on 
an arc. Because of the rapid temporal variability of the burst flux, a 
second stationary scintillator is required to normalize the flux. The 
second configuration considered is an orthogonal pair of stationary 
scintillators of similar dimensions. Both two-detector systems give a 
line position and require a third scintillator to obtain a point position 
for the burst origin. 

A thin flat detector rotating about an axis perpendicular to the 
flux vector will respond to that flux according to a rectified cosine 
function, |cos 0|. If the position of the flux source is not known^the 
function may be fitted to obtain the offset of the source from some 
reference direction. We have performed this fitting for some simulated 
bursts allowing for (a) counting statistics of the quantized flux (b) 
high local and diffuse background (c) intense variation in flux with 
time over very short periods. 


Photon count rates in the detectors were estimated from published 
burst spectra 0 and from natural background and used to estimate the 
location position resolution for various burst strengths* This 
lution was compared with that obtained by fixed orthogonal detectors of 


similar size. 



ROTATING SCINTILLATOR 


For a thin flat scintillator of arbitrary shape, area A cm 2 and 
thickness d cm, the number of counts registered is given by 


N(E,0,t) 



A | cos 0 1 {l-*exp(-ct(E)d sec 0)} 

♦I(E,t)dE d0 dt...(l). 


where 0 is the angle between the flux direction and the detector slab- 
normal, I is the flux intensity in photons cm' 2 s" 1 keV^ 1 and a is the 
absorption coefficient for the detector material. The spin rate was 
chosen to be 2 revolutions s~ T with counts stored in 1024 bins/it 
radians. Dwell time per angular bin is 0.244 ms. Counts per angular 
bin are summed throughout the burst, e-g. for a burst length of 20s, the 
count rate in a particular angular segment is sampled 80 times. Each 
sum of samples for each angular position is stored in one of 1024 bins. 

A similar set of 1024 sums is stored in a corresponding register for 
the stationary detector. Thus the l*-* 1 bin in the register of one 
detector covers exactly the same set of time periods as the iP* 1 bin in 
the second register. 

The burst count SP summed over the burst period for each bin is 
calculated from equation (1). To each of these is added the mean 
background B for the 'same period, which has been previously determined. 
To simulate experimental reality, this set of numbers (SP+B) is 
randomized according to Poisson statistics and is called SP(R). The 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


same process is performed for the stationary detector using a different 
set of random numbers to give ST(R) . 

To remove the effect of temporal variation of burst flux, the 
normalized flux (AA) is calculated from 

AA = -■F.W-.-B-- 
ST(R) - B 

AA still retains the form of (cos 0-[, see figure 1, and is fitted to the 
trial function Kj cos 0) where K is a constant of normalization (figure 2) 
The least squares method then gives the burst direction, Q. This was 
determined ten times using different random number sets and the standard 
deviation of 0 determined for several burst strengths. 

Calculations were performed for a 23 cm x 23 cm x 1 cm Csl(T£) 
crystal at a representative energy of 100 keV. In practice a 1 cm Csl 
slab is rather efficient at this energy and the deviation from the cosine 
function is not great. Edge effects were ignored. 

Several burst strengths were tried and the results of this least- 
squares fit is shown in Table 1. The strongest bursts have a total 
energy of ~ 4 x HT 1 ‘ erg cm" 2 and a total of ~ 400 photons cm"’ 2 burst"” 1 

in the energy interval 40-200 keV are received at the Earth. The burst 
time profile of the April 27, 1972 burst as recorded^ by the Apollo 16. 
shown in figure 3, was used as a model and the count rates scaled up or 
down as required. 

A total estimated background count of 1 photon cm” 2 sec" 1 between 
40 and 200 keV for the local (albedo plus charged-particle-produced) and 
diffuse components was used. This is assumed isotropic and interacts 
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with both front and back surfaces of the detector. This flux is 
typical for a low inclination, low altitude orbit. 



STATIONARY SCINTILLATOR PAIR 


Consider the response of two stationary orthogonal thin flat 
scintillators to a 7-ray flux from a direction perpendicular to the 
line of intersection of their planes. If 0 is the angle between the 
flux vector and the normal to one of the scintillators. 


N. 


0 + tan 1 _A = tan 
Ni 




V B 2 

Ri-Bx 


where R 1 and R 2 are the actual count rati ,s measured for each scintillator, 
and N x 2 and By ? are the respective number of burst and background 
photons respectively for each detector. 

The statistical error in 0, measured over part or all of the burst 
is given by 


cr.(0) = 


n 2 +n 2 
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+N 2 a 2 +N 2 a 2 +n 2 c 2 


j r 2 


R, 




•r 


The position error was calculated for the same burst parameters used 
for the rotating system and is shown in the table. Two cases were 
considered for each burst type; the burst flux vector in the plane of 
one detector, and the flux vector at 45° to both detectors. 




CONCLUSION 


The estimates of angular position errors made here indicate that 
the location uncertanties for gamma- ray bursts cGijld be as low as a 
few minutes of arc for strong bursts and the sizes of scintillators 
considered if systematic errors in measured rates can be kept below 
statistical errors ♦ It is noted that statistical errors fall below 1% 
for the "medium" bursts z\nd these sizes of detector. With respect to 
systematic errors t the rotating scintillator will have an advantage 
over two stationary planar detectors. A direct comparison of two 
detector responses and efficiencies is required for the two fixed 
scintillators whereas for the rotating scintillator a comparison is 
only required to correct for temporal fluctuations in the burst flux. 
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